Polarization of Bloch electrons and Berry phase in the presence of electromagnetic 

fields 
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We consider Bloch electrons in the presence of the uniform electromagnetic field in two- and 
three-dimensions. It is renowned that the quantized Hall effect occurs in such systems. We suppose 
a weak and homogeneous electric field represented by the time-dependent vector potential which is 
changing adiabatically. The adiabatic process can be closed in the parameter space and a Berry 
phase is generated. In the system, one can define the macroscopic electric polarization whose time 
derivative is equivalent to the quantized Hall current and its conductivity is written by the Chern 
number. Then, the polarization is induced perpendicular to the electric field. We show that the 
induced polarization per a cycle in the parameter space is quantized and closely related to the Berry 
phase as well as the Chern number. The process is adiabatic and the system always remains the 
ground state, then, the polarization is quite different from the usual dielectric polarization and has 
some similarity to the spontaneous polarization in the crystalline dielectrics which is also written 
by the Berry phase. We also point out the relation between our results and the adiabatic pumping. 



I. INTRODUCTION 

The presence of the geometrical phase (Berry phase) 
was revealed in the adiabatic process of a quantum me- 
chanical system around a closed loop in a parameter 
space 1^. Thus, it could be regarded as the generalization 
of the Aharonov- Bohm effect |] in a parameter space. 
The Berry phase has been appeared in many contexts, for 
example, a treatment of the Born-Oppenheimer approx- 
imation ^ , fractional statistics , the axial anomaly in 
field theories and so on. 

Recently, the relation between the spontaneous elec- 
tric polarization in crystalline dielectrics and the Berry 
phase was discussed . The authors of Ref . treated 
the system as Bloch electrons with a finite energy gap in 
the absence of electromagnetic fields. They introduce an 
adiabatic change of the potential (Kohn-Sham potential) 
with a slowly varying parameter A. The Hamiltonian is 
compactified in the parameter space, and then, a Berry 
phase is defined in the adiabatic process. The charge 
transfer occurs and the polarization is induced in this 
process with external electric field held to be zero, and it 
is represented by the Berry phase Q . 

The two-dimensional (2D) Bloch electron system with 
an uniform magnetic field has been investigated for a long 
time j^] , and it has been revealed that the integer quan- 
tum Hall effect occurs when the Fermi energy lies in an 
energy gap ||^. The integer part of the quantized Hall 
conductance ne^ /h [n = 0, ±1, ±2 • •■) is represented by 
the Chern number which is a topological number 

defined on the two-torus (the magnetic Brillouin zone). 
The effect is generalized to 3D • Recent arguments 

point out that the band gap could exist for the magnetic 
field around 40 Tesla in organic compounds (TMTSF)2X 
1^,^. Then, to realize the quantum Hall effect of Bloch 
electrons in 3D may be easier than that in 2D. 



In this paper, we consider Bloch electrons in the pres- 
ence of the uniform electromagnetic fields in 2D and 3D. 
A Berry phase is induced by the adiabatic change of the 
time-dependent vector potential. Following Ref. ||^, we 
define the electric macroscopic polarization in the system. 
It is shown that the time derivative of the macroscopic 
electric polarization corresponds to the quantized Hall 
current in the system. The quantized Hall current and 
also the polarization is induced adiabatically. In the adi- 
abatic process, the system always remains to be in the 
ground state, and then, the polarization is quite different 
from the usual dielectric polarization. It has been pointed 
out that the quantized Hall conductivity is written in 
terms of the Berry phase as well as the Chern number 
in the 2D systems [Q. Recently, the authors generalize 
the relation to 3D systems By using the relation 

between the Berry phase and the quantized Hall conduc- 
tance, we can find out that the macroscopic polarization 
is closely related to the Berry phase. The relation be- 
tween the macroscopic polarization and the Berry phase 
seems to be analogous to the spontaneous polarization in 
the crystalline dielectrics 

Recently, the adiabatic pumping is discussed progres- 
sively In pumping, an adiabatic ac perturbation 
yields a dc current, and the charge transfer per the cycle 
is independent on the period of the perturbation. We also 
argue the relation between our results and the adiabatic 
pumping. 

In section II, we consider the 2D system. 3D system is 
discussed in section HI. We set the light velocity c = 1 
throughout the discussion. 
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II. MACROSCOPIC POLARIZATION AND 
BERRY PHASE IN 2D SYSTEM 

In this section, we treat 2D noninteracting electrons in 
a periodic potential in the presence of a uniform mag- 
netic field perpendicular to the plane and a uniform elec- 
tric field in the plane. The electric field is represented by 
a time-dependent vector potential. Then we employ the 
adiabatic approximation assuming that the electric field 
is weak enough. In order to make the paper self-contained 
and fix the notations, we would like to introduce the dis- 
cussions in Ref. first. We review that the existence 
of the Berry phase is demonstrated in this situation and 
that the Chern number is written in terms of the Berry 
phase [p^ . The Chern number is well known as a topo- 
logical number which represents the integer part of the 
quantized Hall conductance in the system |]8|-|l0|]. After 
these reviews, we show that the macroscopic electric po- 
larization, whose time derivative corresponds to the Hall 
current in the system, is written in terms of the Berry 
phase. 

The time-dependent Shrodinger equation we consider 
here is 



(1) 



where 



H{t) 



2m 



(-ihV + eA{t)f + U{v), 



and A(t) is the vector potential for the electromagnetic 
field, therefore, the magnetic field and the electric field 
are given by B = V x A and E — ~dA/dt, respectively. 
We take the Bravais lattice vectors as 



R — ma + Tib, 



(2) 



where m and n are integers, and assume that the poten- 
tial is periodic in both a- and b- directions, i.e. 



U{r) = U{r + a) = U{r + h). 



(3) 



We suppose that the electric field is weak enough. Then 
the system evolves slowly so that it can be described by 
the adiabatic approximation. At any instant, we have an 
eigenvalue equation 



1 



{-ihV + eA{t)f + Uir) 



2m 

£;„(i)*™(i,r) 



(4) 



The system is invariant under translations by a and by 
b. However, the Hamiltonian H(t) is not invariant under 
these transformations. The reason for this is that the vec- 
tor potential A(t) is not constant in space in spite of the 
fact that the magnetic field is uniform. An appropriate 



gauge transformation is required to make the Hamilto- 
nian invariant. It is well known that we can construct 
the magnetic translation operators, which commute with 
Hamiltonian, written as 



Tr = exp 



(5) 



Here we use the symmetric gauge A = (B x r)/2. Now 
we look for eigenstates which simultaneously diagonalize 
Tr and H. However, note that the magnetic translations 
do not commute with each other in general, since 



TaTb = exp[27ri(/)]rbTa, 



(6) 



where (p = {eB/h)\a x b| is the number of magnetic 
flux quanta in a unit cell. When </> is a rational number 
(j) — p/q {p and q are integers with no common factor), 
we have a subset of translations which commute with 
each other. We consider the enlarged unit cell (the mag- 
netic unit cell), in which integral magnetic flux quanta 
go through. For example, if we take the Bravais lattice 
vectors as 



R' ~ ma + nqh, 



(7) 



then p of the magnetic flux quanta penetrates the mag- 
netic unit cell which is enclosed by the vectors a and 
qh. The magnetic translation operators Tr' which corre- 
sponds to the new Bravais lattice vectors commute with 
each other. 

Let be an simultaneous eigenfunction of the oper- 
ators H and Tr/, here a denotes the band index. Then 
the eigenvalues of Ta and T^b are given by 



(8) 



where Tik is the generalized crystal momentum in the 
magnetic Brillouin zone (MBZ). By using the primitive 
vectors of the reciprocal lattice 

2tt 

Ga = — (b X z), 

vq 

Gf, = — (z X a). 



vq ^ z ■ {ax b). 



(9) 



where z is the unit vector perpendicular to the 2D plane, 
k is written as 



k = faGa + -Gb 

q 

< /a < 1, < A < 1. 



(10) 



Therefore, the eigenfunction is labeled by k in addition 
to the band index a. For the magnetic translation invari- 
ance, we can write down the eigenfunction in the Bloch 
form 
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From Eq. 
lation 



4«)(r,i) = 
and Eq.(p 



l"\r,t). (11) 
u'^^r^t) should satisfy the re- 



e u- 



(r + a, t) — exp 



i'^\'r + qh,t) = exp 



i-a-A{t) 



(a) 



U 



ir,t). (12) 



This is the generahzed Bloch theorem in the magnetic 
field 1^. The instantaneous eigenvalue equation for 



u'^\r,t), which is obtained from Eq. (jj) 
is written as 

{-ihV + hk + eA{t)f 



and Eq. (h 



i?k(i) = 



2™ 



C/(r) 



(13) 
(14) 



Now, we have the "Hamiltonian" H]^{t), which contains 
the parameters k and t. Note that the eigenvalue £'^"' (t) 
depends on k continuously. For a fixed band index a, the 
set of values of E^\t) with k varying in the MBZ forms 
the magnetic subband. 

Following Ref. 0] , let us write a time-dependent wave 
function in the adiabatic approximation in terms of 
a snap shot wave function with the dynamical phase 
/q dt'Ek{t') and a phase J^\t) 



*L"^(i)=exp 



E^{t')dt' 



X exp 



(15) 



Substitute the above form into the time-dependent 
Shrodinger equation Eq.(|l|), then, the phase 7^(i) is writ- 
ten as, 



7k (t) = i 



^*(<i>k(i')l^l$k(0)di' 



u^{t')\ — \u^,{t'))dt\ 



(16) 



here we omit the band index a because we will consider 
a single magnetic subband for a while. The "wave func- 
tion" Mk obeys the "Schrodinger equation" (|l3|). 

In order to consider a Berry phase, a Hamiltonian must 
go around a closed loop in the parameter space in the adi- 
abatic process. The Hamiltonian Eq. ( p^ ) as it is does 
not have this property. However, it is possible to com- 
pact ify it as 



(17) 



where the symbol denotes equivalence, since we have 
equivalent crystal momenta and thus give the same wave 
function. The time-dependent part of the vector poten- 
tial is written as —Et. Suppose that the electric field is 



applied along Ga -L b. Then the time-dependence of the 
Hamiltonian //k(i) enters in the form k— (eEvot/hb)Ga- 
Thus we have Uk{t) = Uk-(e£;uot/'ifc)G„, where Uk(t = 
0) = itk, and the period of a transport around a closed 
path in the parameter space is given by T = (hb / cEvq) 
p8[ . The Berry phase is obtained from Eq. (|l^) as 



d_ 

dt' 



"k-f 



d 

' / (Ukl TTT |wk) dfa, 



dt' 



(18) 



where fa and are defined in Eq. (pX]|). Ta{fb) depends 
on /h but is independent on fa- Similarly, if the electric 
field is applied along Gb -L a, we have 



nifa) = t 



/■^ d 
/ {uw\ TTT \uw) dfb 
Jo Ofb 



(19) 



In general, we can define T and Berry phase when the 
electric field is parallel to the reciprocal lattice vector for 
the magnetic unit cell, i.e.. 



E // (mGa + nGb/q), 

where m, n = 0, ±1, ±2 • ••. 

Let us define a vector field in the MBZ by 



A(k) = (uk|Vk|uk), 



(20) 



(21) 



where Vk = 
Then, we have 



dfa. 



dfaGa ■ A(k) 



is the vector operator. 



If dk-A(k), (22) 
'cih) 



where denotes that the path of the line integral is 

taken on the closed loop where fb is fixed. Similarly we 
have 



rb{fa) = i [ dA— ■A(k) 

Jo 1 



dk ■ A(k). (23) 



C{fa 



The vector field A(k) is called Berry connection. This 
is considered as a gauge field induced in the parame- 
ter space, or in the MBZ. Suppose Mk(r) satisfies the 
Schrodinger equation Eq. (13), then an phase trans- 
formed function 



<(r) =Uk(r)e*-^(''), 



(24) 



also satisfies Eq. (p^), where /(k) is the smooth function 
of k and is independent of r. Under the transformation, 
any physical variables do not change, since this transfor- 
mation is global in the coordinate space. From Eq. (I2T 
A(k) is transformed as 



A'(k) = A(k)+zVk/(k). 



(25) 
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Thus, we may regard A(k) as a gauge field induced in 
the parameter space. 

Let us introduce an integral 



^Ch = 



1 

27ri 



MBZ 



Vk X A(k) 



(26) 



Here, [■ ■ ■]z expresses the z-component of the vector. The 
MBZ is topologically a torus rather than a rectangu- 
lar in k-space, because the points k, k+Go and k+Gh/g 
should be identified as equivalent points. Therefore, inte- 
gral coincides with the first Chern number of a principal 
U{1) fiber bundle on the torus whose connection is 
A(k) |p|,p^. Since there are no boundary on the torus, 
the integral (26) becomes finite if and only if A(k) has 
non-trivial topology on the MBZ. 

We apply the Stokes theorem to Eq. ( |26D and derive 
the relation between Berry phases Eqs. (^2|), (^3f) and 
the Chern number Eq. (p™ written as ffil 







dfa 







ajb — -77. 

dfb 



(27) 



Following Ref . Q , one could define the electric contri- 
bution to the macroscopic electric polarization as. 



MBZ (27r) 



^ /ui-\t) r ui-\t)) , (28) 

MBZ ['^^) ^ ' 



where we use Eqs. ( pT[ ) and ([15[), and the summation is 
taken over the bands below the Fermi energy. We rein- 
stalled the band index a. The time derivative of P{t) is 
the electric current of the system. It should be the quan- 
tized Hall current when the Fermi energy lies in the gap. 
Actually, it is calculated as. 



dPjt) 
dt 



MBZ 



+ (4")(t)|r 



d'^k 
(2^ 



d^k 



MBZ 



(2^) 



<'\t)){ui'\t)\r\ui-\t) 



dt 



(29) 



where we use the completeness condition for |Mk(t) > 



Ehfw>(-fw 



and here, the summation is taken over all the bands. We 
can calculate the matrix element for a ^ P 



it) 



it) 



f (t) [Hit),r] ui"\t) 



= —ih- 



(uf (t)|v(t)|4")(t) 



{a) 



e: 



if)) 



E, 



(a) 



it) 



dk 



(30) 



where v(t) ~ {—ihV + eA{t))/m is the velocity operator. 
Then, we can find that 



dP 

'dt 



u 



MBZ 



d'^k ^ 

a<Q:F _ 



it) 



dt 



du'C\t) 
dk 



it) 



dk 



dut\t) 
dt 



(31) 



We have an relation \u\^{t) >~ |wk-eEt/'i >i tbas 
duk{t) eE dui^{t) 



dt 



dk 



Therefore, the time derivative of the macroscopic 
ization is written as 



dP{t) 



polar- 




where denotes the Chern number for the a-th band. 
Eq. ( p^ ) is just the quantized Hall current in the system 
when the Fermi energy is located between the energy gap. 

Then, the polarization induced by the electric field 
along Gq -L b per the period T = (hb/eEvo) (See, the 
paragraph before Eq. (Esl)) is written as 



4 



Ch 



X 



1 

2^ 



ditto 

"/a 







^ ..ditto 
dfb 



, (33) 



where r^"^ is the Berry phase for the a-th band and we 
use Eq. (p7|). The result shows that the electric dipole 
moment per the magnetic unit cell is quantized as an 
integer multiple of — egb. 

Similarly, when E/ /G^ ± a, the period T = ha/eEqvQ 
and 



Q(<;Q:f 



(a) 
Ch ■ 



(34) 



It is also written by the Berry phase by using Eq. (27) as 
well as the Chern number. The electric dipole moment 
per the magnetic unit cell is integer multiple of ea. 

One could point out the specific properties of the po- 
larization in our system as follows; 

(a) introduced by the adiabatic process 

(b) perpendicular to E 

(c) written by the Berry phase 

(d) independent of T, i.e. independent of |E| 

(e) the induced dipole moment per the magnetic unit 
cell per the period T is quantized 



From 



and 



one can see that the polarization in 
our system is quite different in the usual dielectric po- 
larization in the insulators. In the adiabatic process, the 
system always remains in the ground state. In insula- 
tors, the longitudinal conductivity and the translational 
conductivity is zero, and the charge transport would not 
occur adiabatically. Then, the finite polarization is ob- 
tained beyond the adiabatic approximation and the di- 
rection of the polarization is along the external electric 
field E. In our system, we have a finite transverse conduc- 
tivity because of the presence of the magnetic field. The 
transverse charge transport occurs even in the adiabatic 
process and the polarization is induced perpendicular to 
E. 

The properties (a) and (c) suggests that the polariza- 
tion in our system has some similarity to the spontaneous 
macroscopic polarization in the crystalline dielectrics. It 
has been pointed out by King-Smith and Vanderbilt, and 
Resta Q that the spontaneous polarization is induced by 
an adiabatic change of the potential (Kohn-Sham poten- 
tial) with a slowly varying parameter A. The process is 



defined on the closed loop in the parameter space and 
the Berry phase is generated, and the spontaneous po- 
larization is written by the Berry phase. A related result 
is obtained in a different context |l9[ . 

We see the relation between our discussion and the 
adiabatic charge pumping originally argued by Thouless 
The adiabatic pumping is the charge transport 
whose specific properties are: (i) the adiabatic ac pertur- 
bation yields the dc current, (ii) the charge transfer does 
not depend on the period of the perturbation . 

In our system, the perturbation (the electric field) it- 
self is not ac, but the Hamiltonian Hk{t) is compactified 
in the parameter space and the period T is introduced, 
and the dc Hall current Eq. (|32|) flows. 
The statements 



and KeJ imply the fact that the 
charge transfer per T across the boundary of the mag- 
netic unit cell does not depend on T, which corresponds 
to the property (ii) of the adiabatic pumping, and is 
quantized. Actually, one can see it directly. For E ± b, 
the Hall current dP^/dt flows parallel to b and across 
the boundary a-axis. The normal vector of the boundary 
is Gq/|Gq|. Then, the charge transfer is 



AQ 



dt / dfa 



G„ ap. 





IP 



APb 



dt 



For E ± a, the charge transfer is 



AQ = e 



N, 



(a) 
Ch ■ 



(35) 



(36) 



Similar results are obtained by Thouless for the commen- 
surate periodic ac perturbation [ p^ . The commensura- 
bility seems to correspond to the direction of E in our 
discussion (See. Eq. (pO|)), which is essential to define T 
and AQ. 



III. MACROSCOPIC POLARIZATION AND 
BERRY PHASE IN 3D SYSTEM 



The aim of this section is to argue the 3D generaliza- 
tion of the discussions in the previous section. We define 
the macroscopic polarization. One can see that its time 
derivative is equivalent to the quantized Hall current in 
3D. The Hall conductivity is represented by the Chern 
number and quantized when the Fermi energy lies on the 
band gap pd| ]. Recent arguments point out that such a 
gap could exist for the magnetic field around 40 Tesla in 
organic compounds (TMTSF)2X Then, it may 



5 



be easier to realize the quantum Hall effect of Bloch elec- 
trons in 3D rather than 2D. Our discussion written be- 
low could be applicable to such compounds by using the 
tight-binding scheme. 

The relation between the Berry phase and the Chern 
number is also driven in the system ||l5[| . Then, we would 
find out the relation between the macroscopic polariza- 
tion and the Berry phase in 3D system as well as 2D 
system. We also point out the relation between our re- 
sults and the adiabatic pumping fl^Jl^ . 

Let us consider the 3D Bloch electrons in the presence 
of uniform electromagnetic fields. We have a periodic 
potential U{r) = U{r + la + mh + nc) {l,m,n; integers). 
The electric field we consider here is weak and described 
by an adiabatically changing vector potential. We use 
the adiabatic approximation and first, we consider the 
eigenstates of the Hamiltonian H{t) at fixed t. 

We have a magnetic field written as, 



B = B,a + Bbh + B,c, 



(37) 



where a, b and c are the primitive vectors for the Bravais 
lattice. Let Gq = (27r/wo)(b x c), Gb = (27r/'yo)(c x a) 
and Gc = (27r/wo)(a x b) stand for the fundamental re- 
ciprocal lattice vectors, where = a • (b x c). Ba, Bb 
and Be are written as. 









1 




Ba = 


^B. 


Ga = 




-<Pa, 




27r 




vq 


e 








1 


''^ 


Bb = 


±B 


Gb = 








27r 




vo 


e 








I 




Bc = 


±B 


Gc = 




-<Pc, 




2n 






e 



(38) 



where (pa, (pb and (j)c is a number of the unit flux through 
the plane b x c, c x a and a x b, respectively. We in- 
troduce a "rational magnetic fields", where 4>a = Pa/<la, 
4>b = Pb/<lb and 4>c = Pcl<lc are rational values. Here 
Pi, and qi (i = a,b,c) are the integers without common 
factor. Let q stands for the least multiple factor for the 
three integers qa, qb and qc- We can write (pa = Na/q, 
(t>b — Nb/q and (pc — Nc/q {Nf, integer, i = a,b,c). Let 
p stands for the largest common factor for three integers 
Na, Nb and Nc- The magnetic field is written as 

B=l^^c', 



vo eq 

c' = (?^aa + Ubh + iicC) 



(39) 



By definition, Ua, Ub and n,, are integers with no common 
factor. It means that there are no vectors on the Bravais 
lattice which is the submultiple of c', and it was shown 
by Ref. |l^ that we can find vectors a' and b' such that 
a', b' and c' are a new set of primitive vectors as follows. 
Let r stands for the greatest common factor of ric and 
Ua- Therefore, there is no common factor of r and Ub. 



Then we can choose four integers Sa, Sb, Sc and s^, which 
satisfy relations 

Sana + ScTlc = T, SbUb + SrT = 1, 

and we may choose 

a = s^a — SqC, 

, , , SbUc SbUa 

a = s,.b c a. 

r r 

The Hamiltonian has the magnetic translation symme- 
try. Actually, we can choose the magnetic unit cell with 
the primitive vectors a', gb' and c', and the magnetic 
translation operators for the Bravais lattice 

R' = la! + mqh' + nc' {I, m, n; integer), 

is written by the 3D generalization of Eq. (^ when we 
take the symmetric gauge. Three operators T^i, Tq\3i and 
Tc' commute with the Hamiltonian at any instant, and 
also with each other. The wavevector in the magnetic 
Brillouin zone (MBZ) is written as 



fa'Ga' 



q 



Gb'+fc'Ge', 0<fa'Jb'Jc'<l. (40) 



Then, the eigenstates of the Hamiltonian at any in- 
stant is the generalized Bloch states. The Bloch wave 
function with a band index a is written as $^ (r, t) = 

e'^^''^u^\r,t), where u^{r,t) satisfies 



i-iV + -k + eA{t)j^ + U{r) 



1 

2rn 



A{t) = -Et+ -B X r, 



(41) 



which is the 3D generalization of Eqs. ( |T^ ) and ([1^). The 
function u'^\r,t) has relations in the symmetric gauge: 



'(r -f a', i) = exp ■ A{t) u'^^> (r, t), 



<^(r + gb',i) = exp [i^qh' ■ A{t)\ <^r,i), 
4"'(r + c', t) = exp [i Jc' • A(<)] u'^^ (r, t). 



As well as 2D, we have a relation 



(42) 



(43) 



According to Berry ||l|, in the adiabatic approach, 
the solution for the time-dependent Shrodinger equation 
ihd'i/ (t) / dt = H{t)'^{t) can be written with a phase 
ll"\t) as 
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exp 
exp 



-i f dt'E'^\t') 
Jo 



„ik-r„,(a) 



(44) 



We will discuss the presence of the Berry phase later, 
and first, we show the time derivative of the macroscopic 
polarization is equivalent to the quantized Hall current 
in 3D ijTll as follows. The polarization is defined by the 
3D generalization of Eq. (p8|). The definition agrees with 
that in Ref . . Generalizing the derivation of Eq. (|3^ 
to 3D case, we obtain the equation 



dP 

'dt 



-E X D, 



(45) 



where we introduce a "gauge field" defined as A(")(k) 



/, which is the 3D extension of Eq. (25). 
What we have to do here is to investigate the form for 
the three dimensional vector D. We expand D by the 
fundamental reciprocal vectors Gq/, Gb' and Gc'. The 
coefficient of the vector Gr' is written as 



1 

2^ 



(c'.D) 



E 



dh 



Vk X A(")(k) 



c 

1^ 



(46) 



Here, we use the relation Smbz'^^^ ~ /o^ '^/sC^c • 
c'/|c'|) fg^j^jd'^k , where ^'(/a) denotes the plane be- 
tween two fundamental reciprocal vectors G^' and G;,' 
(See, Eq. ^)). The r.h.s. of Eq. {Mj is essen- 
tially equivalent to the Chern number We can 
show in the same manner that ta' = a'-D/27r and 
tb' = q{h' ■ D)/27r are also equivalent to the Chern num- 
ber. Because of the topological nature, these three num- 
bers ta' hi and tc' become integer whenever the Fermi 
energy lies in the energy Gap, even if the magnetic field 
is not rational. By choosing a different MBZ, in which 
the roles of a' and b' are exchanged, we can find that 
b' • D/27r is itself an integer. Therefore, D is a vector 
on the reciprocal lattice generated by Gq' , Gf and Gc' 
written as 



D 



-{ta'Ga' + tb'Gh' + tc'Gc'), 



I.e. 



dPit) 



-E X G. 



(47) 



(48) 



dt 2TTh 

Eq. (^) and (||) suggest that dP/dt is actually the 
quantized Hall current in 3D . 

Then, we will derive the relation between the Berry 
phase and the Chern number in 3D. This is the 3D gen- 
eralization of the discussion in Ref. p4|. Recently, the 
calculation is noted by the authors |15|. 



By substituting Eq. 
Schrodinger equation 



i) into the time-dependent 



J^{t) is obtained as 

li"\t)=^Jyt'(^ui"\t') 



dt' 



'■W-eEt'/h 



dt' 



(49) 



In order to consider the Berry phase, a Hamiltonian 
must go around a closed loop in a parameter space in 
adiabatic process. The Hamiltonian H^it) as it is not 
have this property. However, it is possible to compactify 
it into the magnetic Brillouin zone as -ff/1+1,/2,/3 ~ 
Hfuh+ijAt) ~ ^^/i,/2,/3+i(0 ^ Hf,j^^{t), where /i, 
/2 and /a parameterizes k as Eq. (^) . Repeating the 
discussions in 2D system and by using Eq. (^l|), Berry 
phases in the case that E/ /Ga', E//Gb' and E/ /Gc', are 
written as 



ri?^(/fc',/c') 



ri?\ f c', fa') =^ 



If dk-A(")(k) 

I f rf/„'Ga' • A(")(k), 

Jo 

dk- A(")(k) 

c'Ja') 

.£^G,'.A(")(k), 



(50) 



(51) 



and 



^i"\fa',fb') 



If dk-A(")(k) 

i / d/c'Gc'-A(")(k), 
Jo 



(52) 



respectively. Here, ^^^p f .y (hj = a', b' , d) denotes that 
the path of the line integral is taken on a loop where pa- 
rameters fi and fj are fixed. In general, we can define 
the Berry phase when the electric field is parallel to the 
reciprocal lattice vector for the magnetic unit cell, i.e., 



E // (ZG. + mGb/g + nGc), 



(53) 



where l,m,n — 0, ±1, ±2 • ••. 

We apply the Stokes theorem in Eq. ( p6[ ) and one of 
the three integers tc' , which is equivalent to the Chern 
number, can also be written in terms of Berry phase as 
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a<aF 
Jo "Jb' 



(54) 



and also ta' , tb' are written as 



= ^ E f df-' \ f dfb'^T^:\fa'.h')- 
Jo OJd 

a.<a.Y 

f./.^r(f)(/.,/.) . 



(55) 



(56) 



Then, we calculate the macroscopic polarization change 
per a cycle. As well as 2D, we can define the cy- 
cle when E is parallel to the reciprocal lattice vector 
IGa + rnGb + nGc- For example, the period is written 
as (See, the paragraph before Eq. (p^). 



T 





eE 


h 


|c' xa'l 


qvo 


eE 


h 


a'xb' 



(57) 



Then, the polarization change per T is given; 



-^{qh'tc' - c'qtb') for E ± b'c' - plane, 

' {c'ta' - a'ic' ) for E ± c'a' - plane, (58) 

■^{a'qtb' - qh'ta') for E ± a'b' - plane. 



AP 



For Eqs. (|5J),(^5|) and (|56|), the polarization is written 
in terms of the Berry phase. The dipole moment per the 
magnetic unit cell i.e. gwoAP is quantized in each cases. 

Therefore, one can see that the polarization in 3D also 
has the properties |(aj, |(b)| , |(c)| |(d)| , and (c) written in the 
previous section. Then, the polarization is quite different 
from the usual dielectric polarization and has similarity 
with the spontaneous polarization in dielectric crystalline 
As a property that it is absent in 2D, one can point 
out that the polarization changes its direction discretely 
in the 2D plane perpendicular to E, when one changes 
the location of the Fermi level. The fact reflects the fea- 
ture of the 3D quantized Hall effect fill] . 

As well as 2D, the quantization of the dipole moment 
per the magnetic unit cell per T implies the quantized 
charge transfer across the boundary of the magnetic unit 
cell per T. As the simple extension of Eq. (|35| ) to 3D, 
one can obtain the charge transfer 

-e{qtb' — tc') for E ± b'c' — plane, 
AQ = <( -e{tc' - ta') for E _L c'a' - plane, (59) 
'e{ta' — qtb') for E _L a'b' — plane. 



We should note the fact that the charge transfer is caused 
by the dc Hall current and the results does not depend 
on T. Then, the results is similar to the adiabatic pump- 
ing with the commensurate periodic ac perturbation [ p^ . 
The comniensurability seems to correspond to the direc- 
tion of E in our discussion (See. Eq. (|53|)), which is 
essential to define T and AQ. 



IV. SUMMARY 

In this paper, we have considered Bloch electrons in 
the presence of uniform electromagnetic fields in 2D and 
3D. The Berry phase has been induced by the adiabatic 
change of the time-dependent vector potential. Following 
Ref. [H, the electric macroscopic polarization has been 
defined in the system. It has been shown that the time 
derivative of the macroscopic electric polarization corre- 
sponds to the quantized Hall current in the system whose 
conductivity is represented by the Chern number. Re- 
cently, it was shown that the Hall conductivity is written 
by the Berry phase as well as the Chern number in the 2D 
systems , and also in 3D systems . Then, we have 
found out that the macroscopic polarization is closely re- 
lated to the Berry phase. The quantized Hall current and 
also the polarization has been induced adiabatically. In 
the adiabatic process, the system always remains in the 
ground state, and then, the polarization in our system 
is quite different from the usual dielectric polarization. 
The relation between the macroscopic polarization and 
the Berry phase is analogous to that between the spon- 
taneous polarization in the crystalline dielectrics and the 
Berry phase We have also argued the relation be- 
tween our results and the adiabatic pumping, which is 
discussed progressively fl^Jl^ . 

We have discussed the analogy among the polariza- 
tion in the quantum Hall system, the spontaneous polar- 
ization in dielectric crystalline and the adiabatic pump- 
ing. As an essential point, the analogy comes from the 
fact that these effects are caused by the closed adiabatic 
change in the Bloch electron systems with the finite en- 
ergy gap. 
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